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Abstract — We consider the multiple unicast problem under 
network coding over directed acyclic networks when there are 
two source-terminal pairs, si ti and S2 t2. Current charac- 
terizations of the multiple unicast capacity region in this setting 
have a large number of inequalities, which makes them hard to 
explicitly evaluate. In this work we consider a slightly different 
problem. We assume that we only know certain minimum cut 
values for the network, e.g., mincut(S'i, T,), where Si C {si,S2} 
and Tj C {ti,t2} for different subsets Si and Tj. Based on 
these values, we propose an achievable rate region for this 
problem based on linear codes. Towards this end, we begin 
by defining a base region where both sources are multicast 
to both the terminals. Following this we enlarge the region by 
appropriately encoding the information at the source nodes, such 
that terminal ti is only guaranteed to decode information from 
the intended source Si, while decoding a linear function of the 
other source. The rate region takes different forms depending 
upon the relationship of the different cut values in the network. 

I. Introduction 

The problem of characterizing the utility of network coding 
for multiple unicasts is an intriguing one. In the multiple 
unicast problem there is a set of source-terminal pairs in 
a network that wish to communicate messages. This is in 
contrast to the multicast problem where each terminal requests 
exactly the same set of messages from the source nodes. The 
multicast problem under network coding is very well under- 
stood. In particular, several papers |[T|||2l|l3l discuss the exact 
capacity region and network code construction algorithms for 
this problem. 

However, the multiple unicast problem is not that well un- 
derstood. A significant amount of previous work has attempted 
to find inner and outer bounds on the capacity region for a 
given instance of a network. In an information theoretic 
characterization for directed acyclic networks is provided. 
However, explicit evaluation of the region is computationally 
intractable for even small networks due to the large number of 
constraints. The authors in Q propose an outer bound on the 
capacity region. Price et al. [jS] provide an outer bound on the 
capacity region in a two unicast session network, and provided 
a network structure in which their outer bound is the exact 
capacity region. The work of |i7J forms a linear optimization 
to characterize an achievable rate region by packing butterfly 
structures in the original graph. This approach is limited since 
only the XOR operation is allowed in each butterfly structure. 
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In this work we propose an achievable region for the two- 
unicast problem using linear network codes. Our setup is 
somewhat different from the above-mentioned works in that 
we consider directed acyclic networks with unit capacity edges 
and assume that we only know certain minimum cut values 
for the network, e.g., mincut(S'i, Tj), where Si C {si, S2} and 
Tj C {ii, <2} for different subsets Si and Tj. This is related to 
the work of Wang and Shroff (Sj (see also |j9j|) for two-unicast 
that presented a necessary and sufficient condition on the net- 
work structure for the existence of a network coding solution 
that supports unit rate transmission for each Si — ti pair. In this 
work we consider general rates. Reference ifTOl is related in the 
sense that they give an achievable rate region for this problem 
based on the number of edge disjoint paths for Si — ti pair 
In our work we propose a new achievable rate region given 
additional information about the network resources. The work 
of ifTTl considered the three unicast session problem in which 
each source is transmitting at unit rate. Finally, reference lfT2l 
applies the technique of interference alignment in the case of 
three unicast sessions and shows that communication at half 
the mincut of each source-terminal pair is possible. 

This paper is organized as follows. Section HIl introduces the 
system model under consideration. Section|III]contains the pre- 
cise problem formulation and the derivations of our proposed 
achievable rate region. Section |IV] compares our achievable 
region to existing literature. Due to space limitations, some of 
the lemma proofs are not given and can be found in |fT3l . 

II. System Model 

We consider a network represented by a directed acyclic 
graph G = (V, E). There is a source set S — {si, S2} G ^ in 
which each source observes a random process with a discrete 
integer entropy, and there is a terminal set T = {^1,^2} £ V in 
which ti needs to uniquely recover the information transmitted 
from Si at rate Ri. Each edge e £ E has unit capacity and can 
transmit one symbol from a finite field of size q. If a given edge 
has a higher capacity, it can be divided into multiple parallel 
edges with unit capacity. Without loss of generality (W.l.o.g.), 
we assume that there is no incoming edge into source Si, and 
no outgoing edge from terminal ti. By Menger's theorem, the 
minimum cut between sets Sni ^ S and Tjv^ C T is the 
number of edge disjoint paths from Sni to T^^, and will be 
denoted by k^-^^N^ where Ni,N2 C TV ~ {1>2}. For two 
unicast sessions, we define the cut vector as the vector of the 



outvalues fci_i, fc2-2, ^1-2, ^2-1, ^12-1, ^12-2, fci-12, ^2-12 
and fci2-i2. 

The network coding model in this work is based on 
ID. Assume source Si needs to transmit at rate Ri. Then 
the random variable observed at Si is denoted as Xi = 
{Xii, Xi2, ■ ■ ■ ,Xiii.), where each Xy is an element of 
GF{q); the XiS are assumed to be independent. For linear 
network codes, the signal on an edge is a linear combi- 
nation of the signals on the incoming edges on i or a linear 
combination of the source signals at i. Let Y^^ (tail (en) = k 
and head{en) = I) denote the signal on edge e„ G E. Then, 

Ye„ = fm,nYe^ if k E V \ {si, S2}, and 

{cin I head{ejn)—k} 

Ri 

Ye^ = aij,nXij if Xi is observed at k. 

i=i 

The local coding vectors aij^n and fm.n are also chosen from 
GF{q). We can also express Ye^ as , Yg^ — J^fli ctj.nXij + 
^f=i P 3,^X2 j- Then the global coding vector of Yg^ is 

= [«!,„, a2,n, ■• ■ ,a_Ri,ra,/3l,n,/32,Ti, • • • ,/3fl2,n]- We 

are free to choose an appropriate value of the field size q. 

In this work, we present an achievable rate region given a 
subset of the cut values in the cut vector; namely, fci^i, ^2-2, 
fci-2, ^2-1, fci2-i, fci2-2- W.l.o.g, we assume there are fci_ij 
outgoing edges from Si and kij^i incoming edges to t^. If this 
is not the case one can always introduce an artificial source 
(terminal) node connected to the original source (terminal) 
node by ki^ij (kij^i) edges. It can be seen that the new 
network has the same cut vector as the original network. 

III. Achievable rate region for a given cut vector 

First, suppose that only ti is interested in recovering the 
random variables Xi and X2 which are observed at si and 
S2 respectively. Denote the rate from si to ti and S2 to ti as 
Rii and Ri2. Then the capacity region Ct^, that is achieved 
by routing will be 

Rii < fci-i, 
R12 < k2-i, 

Rii + R12 < fcl2-l- 

The capacity region Ct2 for t2 can be drawn in a similar 
manner This is shown in Fig. |l(a)| We also find the boundary 
points a,b,c,d such that their coordinates are a = (fci2-i — 

^2-1,^2-1),^ = (fcl-2,fcl2-2 - fcl-2),C = (fci_i,fci2-l - 

= (fci2-2 — ^2-2,^2-2)- A simple achievable rate 
region for our problem can be arrived at by multicasting both 
sources Xi and X2 to both the terminals ti and ^2- 

Theorem 3.1: Rate pairs (_Ri, R2) belonging to the follow- 
ing set B can be achieved for two unicast sessions. 

B = {Rl < min(fci_2, fci-i), 
i?2 < min(fc2_i, ^2-2), 
i?i + i?2 < min(fci2-i, fci2-2)}- 

Proof: We multicast both the sources to each terminal. 
This can be done using the multi-source multi-sink multicast 
result (Thm. 8) in 12. ■ 



R2 




Fig. 1 . (a) An example of a capacity region, (b) Base region for the example. 

Subsequently we will refer to region B achieved by multi- 
cast as the base rate region (the grey region in Fig. |l(b)| i. 

We now move on to precisely formulating the problem. Let 
Zi denote the received vector at ti, Xi denote the transmitted 
vector at Si, and Hij denote the transfer function from Sj 
to ti. Let Mi denote the encoding matrix at Si, i.e.. Mi is 
the transformation from Xi to the transmitted symbols on 
the outgoing edges from s^. In our formulation, we will let 
the length of Xi to be fcj_, (i.e., the maximum possible). 
For transmission at rates i?i and R2, we introduce precoding 
matrices Vi,i = 1, 2 of dimension Ri x ki-i, so that the overall 
system of equations is as follows. 

Zi = HuMiViXi + H12M2V2X2, 
Z2 = H21M1V1X1 + H22M2V2X2. 

We say that ti can receive at rate Ri from Si if it can decode 
ViXi perfectly. The row dimension of the K/s can be adjusted 
to obtain different rate vectors. For (i?i,i?2) £ B, it can be 
shown that there exist local coding vectors over a large enough 
field such that the ranks of the different matrices in the first 
column of Table H] are given by the corresponding entries in 
the third column, which correspond to the maximum possible. 
Furthermore, by the multi-source multi-sink multicast result, 
these matrices are such that [HuMi H12M2] is a full rank 
matrix of dimension fci2-i x (_Ri+i?2), and [H2l^^l H22M2] 
is a full rank matrix of dimension fci2_2 x (_Ri + R2). In 
Table H] for instance since the minimum cut between si 
and ti is fci_i, we know that the maximum rank of Hu 
is Using the formalism of |2j, we can conclude that 

there is a square submatrix of Hu of dimension fci_i x fci_i 
whose determinant is not identically zero. Such appropriate 
submatrices can be found for each of the matrices in the first 
column of Table I] This in turn implies that their product is 
not identically zero and therefore using the Schwartz-Zippel 
lemma, we can conclude that there exists an assignment of 
local coding vectors so that the rank of all the matrices is 
simultaneously the maximum possible. For the rest of the 
paper, we assume that such a choice of local coding vectors has 
been made. Our arguments will revolve around appropriately 
modifying the source encoding matrices Mi and A/2. 

Note that there are two boundary points of the base region 
(the two boundary points may overlap). At point Qi, we 
denote the achievable rate pair by (i?*, i?2) where 

i?^ = min(fci_2, and 

i?2 = min(min(fc2_i,fc2-2),min(fci2_-i,fci2-2) - R*)- 



TABLE I 
DIMENSION OF MATRICES 



matrix 


dimension 


rank 




fcl2-l X fcl-12 


ki-i 




fel2-l X fc2-12 


^2-1 


[Hii H12] 


kl2-l X (fci_i2 + fe2-12) 


^12-1 


Ml 


fcl_12 X Rl 


Bi 


H21 


fel2-2 X fel-12 


ki-2 


H22 


fel2_2 X fc2-12 


k2-2 


[H21 H22] 


ki2-2 X (fci_i2 + fe2-12) 


ki2-2 


Ah 


fe2-12 X ri.2 


R2 



At point Q2, we denote the achievable rate pair by {R\* ^R^*) 
where 

i?2* = min(fc2_i,fc2_2), and 

Rl* = miii(min(fci_2, fci-i),min(fci2_i, fci2-2) - -R2*) 



In Fig. |l(a)| these boundary points are Qi = 6 and Q2 = e. 

In what follows, we will present our arguments towards 
increasing the value of i?i to be larger than Rl (these 
arguments can be symmetrically applied for increasing R2 as 
well). For this purpose, we will start with the point Qi and 
attempt to achieve points that are near it but do not belong to 
B. At Qi, if Rl = fci_i, then we cannot increase Ri due to the 
cut constraints. Hence, we assume Rl ~ fci_2- Furthermore, 
since fc2-2 > ^12-2 - fci-2 > min(fci2-i, fci2-2) - fci-2, 
i?2 = min(min(fc2-i, fc2-2), min(fci2-i, fci2-2) - -Rl) = 
min(fc2-i, min(fci2_i, fci2-2) - ^1-2)- 

In this paper we refer to fci_2 + k2-i as a measure of the 
interference in the network and in the subsequent discussion 
present achievable regions based on its value. We emphasize 
though that this is nomenclature used for ease of presentation. 
Indeed a high value of ki-2 does not necessarily imply that 
there is a lot of interference at t2, since the network code itself 
dictates the amount of interference seen by t2- The following 
lemma will be used extensively. 

Lemma 3.2: Consider a system of equations Z = HiXi + 
H2X2, where Xi is a vector of length li and X2 is a vector 
of length I2 and Z G span{[Hi H2]j^ The matrix Hi has 
dimension zt x li, and rank /i — cr, where < <t < li. 
The matrix H2 is full rank and has dimension zt x I2 
where zt > {h + h — cr). Furthermore, the column spans 
of Hi and H2 intersect only in the all-zeros vectors, i.e. 
span{Hi) n span{H2) = {0}. Then there exists a unique 
solution for X2. 



A. Low Interference Case 
This is the case when fci- 



+ fc2~i < min(fci2-i, fci2-2) 



At Qi, from the assumption, it follows that Rl = fci_2, -R2 ^ 
min(fc2-i, min(fci2-i, fci2_2) - fci-2) = h-i- An example 
is shown in Fig. |2(a)| Furthermore, Qi = Q2 = e. 

'Throughout the paper, span{A) refers to the column span of A. 



Our solution strategy is to consider the encoding matrices 
Ml and M2 at the point Qi, and to introduce a new encoding 
matrix at si, denoted M[ (with Rl + 6 columns) such that 
span{HiiM[)f\span{H 12) = {0}. As shown below, this will 
allow ti to decode from si at rate Rl+5 and t2 to decode from 
S2 at rate i?2. After the modification, each ti is guaranteed to 
decode at the appropriate rate from si. A similar argument can 
then be applied for i?2 to arrive at the achievable rate region 
in this case. 





Fig. 2. (a) The capacity regions Cti and Ctj for an example of low 
interference case, (b) The achievable rate region for low interference case. 
For each point in the shaded grey area, both terminals can recover both the 
sources. In the hatched grey area, for a given rate point, its x-coordinate is 
the rate for si — ti and its y-coordinate is the rate for S2 —t2; the terminals 
are not guaranteed to decode both sources in this region. 

At the point Qi, the rates are Rl = fci_2, R2 = ^2-1- Since 
both terminals can decode both sources, it holds that 

rank{HiiMi) ~ ki^2,'i'o.nk{Hi2M2) = fc2-i, and 
span{H,iMi) n span{Hi2M2) = {0} for i ^ 1,2. 

By analyzing the properties of the above matrices, we have 
Theorem 13.41 Before we state the theorem, we first give the 
following lemma which will be used in proving Theorem 13.41 

Lemma 3.3: Rate Increase Lemma. In the base region, 
denote the achievable rates at Qi as Rl and R2, and 
the corresponding encoding matrices as Mi and M2. Let 
rank{[Hii H12M2]) =^ r > Rl + RX. There exist a series 
of full rank matrices M^^"'' = [m["'' Mi] of dimension 
fcl-12 X (n + Rl) such that rank{[HiiM['"^ H12M2]) = 
Rl+R^+ n, < 71 < (r - i?* - R*). 

Theorem 3.4: Given a cut vector, if fci_2 + fe-i < 
min(A:i2-i, ^12-2), then the rate pair in the following region 
can be achieved. 
Region 1: 



^1 < ki2- 
R2 < ki2- 



^2-1, 
fcl-2, 



which is shown in Fig. |2(b)| 

Proof: In this case, Rl = ki-2 and i?2 = ^2-1 is the 
boundary point Qi = Q2. We will try to find full rank matrix 
Ml of dimension A:i_i2 x (A:i2_i — ^2-1) and full rank matrix 
M2 of dimension fc2-i2 x (^12-2 
of equations can be written as 

Zi = HiiM[V;Xi H 
Z2 = H2iM[ViXi H 



ki-2) such that the system 



Hi2M'2ViX2, 
H22M'^V^X2, 

and VlXi can be decoded at ti, V2X2 can be decoded at t2. 



First, note that rank{Hi2M2) = rank{Hi2), which 
implies that span{Hi2) = span{Hi2M2)- Therefore 
rank{[Hii H12]) = rank{[Hii H12 H12M2] = 
rank{[Hii Hi2M2])- Together, this implies that 
rank{[Hii Hi2M2\) — fci2_i. Using the Rate Increase 
Lemma, we can find the matrix AI[ such that the 
following conditions are satisfied: (i) M[ is a full 
rank matrix of dimension fci_i2 x (fci2-i — fc2-i), 
(ii) rank{HiiM[) = fci2-i — ^2-1 and (iii) 
span{HiiM[) D span{Hi2) = {0}. (i) is from the Rate 
Increase Lemma, (ii) and (iii) hold because of the following 
argument. From Rate Increase Lemma and the fact that 
rank{Hi2M2) = rank{Hi2) = k2-i, we will have 

fci2-i = rank{HiiM[ H12M2) 

= rank{HiiM[) + rank{Hi2M2) 

— rank{span{Hii M[) fl span{Hi2M2)) 

< rank{HuM[) + rank{Hi2M2) 

< rank{M[) + rank{Hi2) 

= fci2-i - ^2-1 + ^2-1 = fci2-i- 

Then all the inequalities become equalities, (ii) and (iii) are 
satisfied. Likewise, can be found with similar conditions. 

Next, since span{HiiM[) D span{Hi2) = {0} and 
span{Hi2M2) C span{Hi2), we will have span{HiiM[) n 
span{Hi2M2) = {0}. By Lemma 13.21 and the above three 
conditions, ti can decode V{Xi at rate fci2-i — ^2-1, but 
cannot decode V2X2. By a similar argument, t2 can decode 
V2X2 at rate ^12-2 — ^1-2, but cannot decode V{Xi. ■ 



B. High Interference Case 

This is the case when fci_2 + ^2-1 > inin(fci2-i, fci2-2)- 
Recall that we also assume that fci_2 < fci-i. At Qi, 
R* = fci_2, R2 = min(fc2-i, miii(fci2_i, fci2-2) - ^1-2) = 
min(fci2-i, fci2-2) — fci-2- This means that Qi and (52 are 
two separated points. An example is shown in Fig. |l(a)| In 
particular, when Ct^ is contained in Ct2 or vice versa, the 
achievable region is described by this case. 

Our strategy is similar to the one for the previous case, 
but with important differences. We begin with the rate vector 
at point Qi and then attempt to increase However, in 
this particular case we will not be able to increase R2 and 
in fact may need to reduce it. This is because at point Qi, 
we have i?2 = rank{Hi2M2) < ^2-1 = rank{Hi2), i.e., 
the encoding matrix ]\I2 is such that rank{Hi2M2) is strictly 
less than the maximum possible. Therefore, if we augment 
A/2 with additional columns to arrive at Mj, it is not possible 
to assert as before that the span{HiiM[) D span{Hi2M2) = 
{0}. Hence, it may be possible that si cannot be decoded at 
ti, (after augmenting M2 to M2). In this situation, we have 
the following result. 

Theorem 3.5: Given a cut vector, if fci_2 + ^2-1 > 
min(fci2_i, fci2_2) and fci_2 < fci-i, then the rate pair in 
the following region can be achieved. 



Region 2: 

Ri < fci-i, 

i?2 < min(fci2_i, fci2-2) - fci-2, 
i?i + i?2 < rank{[Hn H12M2]). 

Note that in the above characterization, the sum rate constraint 
depends on rank{[Hii H12M2])', we show a lower bound 
on rank{[Hii H12M2]) in IIII-Bll The following lemma that 
discusses situations in which rates can be traded off between 
the two unicast sessions is needed for the proof of Thm. 13.51 

Lemma 3.6: Rate Exchange Lemma. Given that 
rank{[HiiMi H12M2]) = rank{[Hii H12M2]) = r, where 
Ml is a full rank matrix of dimension ki-12 x (r — R2), 
M2 is a full rank matrix of dimension fc2-i2 x i?2- If 
M[ = [ d All ] where o? is a vector of length fci_i2 and 
rank{HiiAI[) = r — R2 + I, then there exists an A/j such 
that span{HiiM[) n span{Hi2M^) = {0} where M.^ is a 
full rank submatrix of M2 of dimension fc2-i2 x {R2 — 1). 
Proof of Theorem 15.51 Given that fci_2 + fe-i > 
min(fci2-i, ^12-2) and fci_2 < ki^i, we will extend the rate 
region from Qi where i?^ = fci_2, i?2 = miii(fci2-i, ^12-2)^ 
fci_2. At Qi, we need to increase Ri while keeping R2 as 
large as possible. By the Rate Increase Lemma, we can achieve 
the rate point R'l = rank{[Hii Hi2M2])-Rl, R'2 = R^- The 
corresponding encoding matrices are M[ and A/2. When we 
want to further increase R'^, we could use the Rate Exchange 
Lemma repeatedly. Hence, when R'l is increased by 5, R'2 
is decreased by 5 where Q < 5 < min(/?2,fci_i — R'l) 
{S < fci-i — R'l comes from the fact that R'l can be increased 
to at most fci-i). Terminal ti can decode messages both from 
si at rate /?" = R'^ + S and S2 at rate R'2 = R'2 — S. Denote 
the new set of encoding matrices as M" and A/j'- 

At t2, because A/2' is a submatrix of A/2, span(i/22^^2') — 
span{H22Ad2) ■ Furthermore, we have span{H2iM'{) C 
span{H2i) = span{H2iMi), since RI ~ fci-2- Hence, 
from the above argument, we will have span{H2iM'{) n 
span{H22Ad'2) = {0} since span{H2iMi)r)span{H22Al2) = 
{0}. Then by Lemma [3721 we can decode at R'2 ^ R'2 — S 
from 32, but not decode any messages from si. ■ 

A similar analysis for Q2 allows us to increase i?2, resulting 
in the following extended region. 

Corollary 3.7: Given a cut vector, if ki-2 + ^2-1 > 
min(/ci2_i, ^12-2) and fc2-i < ^2-2, then the rate pair in 
the following region can be achieved. 
Region 3: 

Ri < min(fci2_i, fci2-2) - ^2-1, 

R2 < fc2-2, 

Ri+R.2< rank{[H2iMi //22]). 

The overall rate region is the convex hull of base region. 
Region 2 and Region 3 which is shown in Fig. |3(a)[ where 
boundary segment d — / is achieved via timesharing. 

We note that the idea of increasing one rate while decreasing 
the other can also be applied to the region obtained in low 
interference case. Since rank{[Hii H12M2]) = fci2-i and 
rank{[H2iMi H22]) = ^12-2, we can obtain the following 
two new regions for low interference case. 



Region 2': 

Ri < ki^i 

i?2 < k2-l 

Ri+R2< ki2-i 



Region 3': 

Ri < fci-2 

i?2 < ^2-2 
i?l + i?2 < A:i2-2 



Finally, the achievable rate region for low interference case 
is the convex hull of the region 1, 2' and 3' shown in Fig. 
|3(b)[ where the boundary segment — / and / — c is achieved 
via timesharing. 

1) Lower bound of rank{[Hix Hi2M2\): Next, we investi- 
gate the lower bound of rank(\Hii Hi2M2\). In the follow- 
ing argument, and i?2 denote the rate at boundary point Qi, 
and Ml and M2 denote the corresponding encoding matrices. 
First note that ranA:([iJii Hi2M2\) > rank^Hu) = /ci-i and 
rank{[Hii H12M2]) > rank{[HiiMi H12M2]) = RI + R^ 
Next we will also find another nontrivial lower bound of 
rank{[Hii H12M2]) by the following lemma. 

Lemma 3.8: Given rank{[Hii ^^12]) = ^12-1, 
rank{Hi2) = ^2-1 and rank{[Hi2M2]) = I, we have 
rank{[Hii H12M2]) > fci2_i - fca-i + I. 

Proof: By the assumed conditions, there are fc2-i 
columns in H12 that are linearly independent, and in Hn, 
we can find a subset of fci2_i — fc2_i columns denoted 
such that span{H[i) D span{Hi2) = {0} and rank{H[i) = 
fci2-i — ^2-1, which further imply that rank{[H[i H12]) = 

fcl2-l- 

Since span{Hi2M2) C span{Hi2) this means 
that span{H[^) n span{Hi2M2) = {0}. Then 
rank{[H[^ H12M2]) = rank{H[^) + rank{Hi2M2) - = 
fci2-i — ^2-1 + I- Hence, rank{[Hii H12M2]) > 
rank{[H[^ H12M2]) = ^12-1 - fca-i + /. ■ 

Together with the two lower bounds above, we have 
rank{[Hii H12M2]) > max(fci_i, fci2-i -fc2-i+^27^i + 
i?2)- A case where max(fci_i, fci2-i — fc2-i +^21 -^i +^2) 
fci2-i — k2-i + i?2 is shown in Fig. |3(a)| 




R2 



(b) 



RI 



Fig. 3. (a) The extended rate region for liigli interference case, (b) The final 
extended rate region for low interference case. For each point in the shaded 
grey area, both terminals can recover both the sources. In the hatched grey 
area, for a given rate point, its z-coordinate is the rate for si — ti and its 
i;-coordinate is the rate for S2 — i2; the terminals are not guaranteed to decode 
both sources in this region. 

IV. Comparison with existing results 

The authors in JS) and JO) explore the case when each source 
transmits one symbol at a time, or equivalently, i?i = i?2 = 1 
in detail, whereas we allow arbitrary rate pairs. Reference 
ifTol , also consider the scenario where the rates are arbitrary. 
Assuming that fc2-2 < ^i-i, the basic region in jTOl is 



Region EF09: 

i?i + 2i?2 < fcl-l 
R2 < k2-2 

They also extend the region using the knowledge of fci_2, fe-i 
and other cut conditions arising from the network topology 
(see section IV of IfTOl ). A comparison between our region 
and theirs indicates that there are example networks where 
there exist rate points that belong to our region but not to 
Region EF09. Conversely, there are instances of networks 
where points that belong to Region EF09, do not fall within 
our region. The work of (TOj can be interpreted in part as 
an interference nulling scheme, and in future work it may 
be possible to incorporate this within our approach. The 
work of [6| considers several different cuts defined in the 
graph and propose an outer bound for the network capacity. 
Moreover, they provide certain network structures where the 
outer bound is tight. Since our work deals with an inner bound, 
it is qualitatively different. Finally Das et al. |[T2| have used 
interference alignment for the case of three unicast sessions, 
and are able to achieve a rate that is half the mincut for each 
unicast session. While this is an interesting result for a harder 
problem, the case of two unicast sessions considered here is 
different since each connection has only one interferer and 
the alignment problem does not exist. Moreover, achieving 
half the mincut for each session can be trivially achieved 
by timesharing in our problem. In that sense a comparison 
between our results and theirs is not possible. 
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